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Abstract 


High  angle  of  attack  flight  control  is  of  utmost  importance  to  military 
aircraft  in  air  combat  maneuvering.  Flight  in  this  regime  has  in  recent 
jears  caused  many  high  performance  aircraft  to  be  lost  due  to  departure  of  the 
aircraft.  The  aerodynamics  in  this  regime  are  highly  nonlinear.  The  problem 
is  compounded  by  the  fact  that  the  aerodynamics  are  not  well  known.  This 
paper  considers  the  use  of  adaptive  control  in  order  to  perform  model  following 
of  an  "ideal"  aircraft  in  the  presence  of  uncertain  aerodynamic  coefficients. 

In  particular,  the  partitioning  approach  of  adaptive  control  is  extended  to 
tiie  implicit  model  following  problem.  This  is  then  used  to  solve  the  problem 
of  high  angle  of  attack  flight  control. 

I . Introduction 

High  angle  of  attack  flight  may  cause  undesirable  aircraft  dynamic 
response.  This  dynamic  response  has  been  the  cause  for  the  loss  of  many 
k fiigh  performance  swept-wing  aircraft  due  to  stall -departure  problems.  Air- 

c^aft  departure  can  be  defined  [1]  as  an  uncommanded  and/or  uncontrollable 
^ dynamic  response  of  the  aircraft.  It  is  manifested  as  either  a divergent 

t ilinq-side  slipping  oscillation  of  large  amplitude,  i.e.,  wing  rock,  or  a 
!a>ge  rapid  yaw  generally  followed  by  a rapid  roll,  i.e.,  nose  slice.  For 
. example,  reference  [1]  documents  a simulation  in  which  an  A-7  airplane  in  a 

turn  at  high  angle  of  attack  flight  exhibited  nose  slice  with  a yaw  rate 
buildup  of  approximately  65  deg/sec.  In  an  actual  flight  this  would  have 
caused  a spin  from  which  the  probability  of  recovery  would  have  been  low. 

'•  When  this  condition  is  encountered,  it  is  unanticipated  and  the  pilot  is  under 

f'hysical  and  mental  stress;  the  aircraft  and  pilot  will  most  likely  be  lost  if 
departure  occurs  during  combat  and  certainly  will  be  lost  under  15,000  ft 
altitudes  regardless  of  spin  recovery  characteristics  [2]. 

The  importance  of  depa>"ture  led  to  a symposium  at  the  Air  Force  Flight 
Ovnamics  Laboratory  in  1971  on  AFFDL  Stall/Post  Stall/Spin  Symposium.  Since 
‘.hen  many  articles  have  appeared  in  the  study  of  departure  [1,2, 3, 4, 5]. 

Nlumerous  studies  have  been  made  to  identify  the  aerodynamic  coefficients  in 
high  angle  of  attack  regime.  Reference  [6]  is  a complete  study  of  this. 


However,  in  order  to  collect  real  time  data  for  proper  identification,  it  is 
necessary  to  subject  the  aircraft  to  divergence.  Wind  tunnel  data,  although 
valuable,  still  have  residual  errors  due  to  differences  in  unsteady  flow 
between  the  wind  tunnel  and  the  actual  flight  condition.  It  is  desirable  to 
accurately  have  knowledge  as  to  the  coefficients  for  control  purposes;  in  fact, 
it  is  necessary. 

Ihis  paper  considers  the  control  of  aircraft  at  high  angles  of  attack  in 
the  presence  of  uncertainty  in  the  aerodynamic  coefficients.  In  particular, 
tiie  nonlinear  equations  of  motion  are  linearized  about  a given  flight  condi- 
tion. The  nonlinear  aerodynamics  are  included  in  the  model.  An  ideal  aircraft 
model  for  the  given  flight  condition  is  developed.  This  ideal  model  varies 
with  charijes  in  flight  condition.  An  adaptive  implicit  model  following  control 
law  is  developed  in  order  to  keep  the  actual  aircraft  close  to  the  ideal 
response.  Uncertainties  in  the  aerodynamic  coefficients  are  eliminated  by 
t*'  adaptive  estimator. 

This  paper  is  divided  into  five  sections.  The  next  section  contains  the 
problem  statement.  .Section  III  gives  the  development  of  the  control  law. 
Section  IV  contains  a discussion  as  to  how  the  ideal  model  for  the  A-7  was 
boson.  A different  aircraft  may  have  a different  "ideal"  model.  Section  V 
c-ntains  simulation  results  for  the  A-7  aircraft  using  the  control  law  in  the 
paper.  Section  VI  yields  the  conclusions. 


1 1 . Problem  Statement 

The  equations  of  motion  of  an  aircraft  linearized  about  a given  angle  of 
attack,  x^-,  Euler  angles  between  gravity  oriented  inertial  axis  and  aircraft 
body  axis,  the  flight  path  angle,  YqI  angular  velocities,  p^,  q^, 

and  r^;  nominal  forward  velocity,  U^;  sideslip  angle,  6^,  and  given  as 


whern 


X = A(uj ,t)x  + 6(^2. T)u 


' |“v  '‘■'"o’  '>■  <'•  “p'*po|  ’ 


' |‘e’  ‘s’  ‘r|' 
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wh.nf?  is  the  perturbed  total  linear  velocity;  m-a^  is  the  perturbed  angle 
of  attack;  p,  q,  and  r are  the  perturbed  angular  velocities;  is  the 

perturbed  sideslip  angle;  *-^p"®po  the  perturbed  pitch  angle  and  '^is  the 
foil  angle.  The  controls  are  assumed  to  be  deflections  in  elevator,  6^,  in 
aileron,  6^,  and  in  rudder,  6n.  The  matrices  A and  B are  given  in  Appendix  A 
as  well  as  the  definitions  of  Uj  and  ^2-  parameter  vector  uj  and  ^2 

contain  the  aerodynamic  coefficients  which  are  assumed  uncertain  except  for  an 
a priori  probability  density  function.  The  coefficients  are  assumed  constant 
over  the  time  interval  that  the  linearization  is  assumed  valid.  The  time 
dependency  of  the  A and  B matrices  depict  the  temporally  changing  lineariza- 
tion. 

It  is  assumed  that  a noisy  measurement  of  the  states  is  available,  i.e., 

= Cx  + V (2) 

where  v is  zero  mean  white  noise  with  covariance  E|v(t)v(T)^j  = Vp6(t-i)  and 
is  defined  in  Appendix  A. 

From  Equation  (A. 2 ),  it  may  be  noted  that  the  longitudinal  and  lateral 
mooes  of  the  aircraft  are  highly  coupled.  Furthermore,  as  may  be  noted  by 
expansion  of  the  equations  of  motion,  there  are  many  destabilizing  terms  in 
the  equations.  Consequently,  it  is  an  extremely  difficult  multivariable  task 
for  Che  pilot  to  prevent  departure  at  high  angles  of  attack.  This  is  especially 
true  in  air  combat  maneuvering  as  the  physical  and  mental  stress  occupies  the 
pilot's  attention.  Thus,  it  is  desirable  to  obtain  a closed  loop  control  law 
which  will  prevent  departure.  One  method  is  reported  in  [5]  where  the  author 
develops  a feedback  control  to  eliminate  perturbations  about  a nominal  trajec- 
tory for  a deterministic  system.  The  approach  in  this  paper  is  to  design  a 
feedback  control  law  based  on  model  following  of  an  ideally  responding  air- 
craft. Ihe  advantage  is  that  the  model  following  control  more  closely  gives  a 
control  law  that  will  yield  a desirable  response.  Thus,  for  example, 
decoupling  of  the  yaw-roll  problem  in  nose  slice  may  be  approximately  obtained 
by  placing  this  feature  into  the  ideal  model.  The  destabilizing  terms  may  be 
compensated  by  use  of  the  model.  The  adaptive  control  law  will  compensate 
for  ur certai nties  in  the  aerodynamic  coefficients  by  real  time  learning.  The 
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development  of  the  ideal  model  will  be  explained  in  Section  IV  with  a specific 
example  given  for  the  A-7  in  a particular  flight  condition.  The  form  of  the 
model  is 


z = \(t)z  (3) 

where  the  time  varying  matrix  corresponds  to  the  ideal  model  changing  due 
to  the  different  flight  conditions.  Since  implicit  model  following  [10]  is  to 
be  accomplished,  the  performance  index  is  taken  to  be 

‘o 

where  is  the  output  vector  (not  to  be  confused  with  the  measurement  vector) 

Vo  = V (5) 

^ where  is  a time  invariant  distribution  matrix  and  is  a general  term 

’ since  x and  z ma:,  not  be  of  the  same  dimensions  and  where  Rp  weights  the 

control  surface  deflections  and  is  a positive  definite  matrix,  Qp  weights 
excursions  from  the  model  response,  and  t^  is  chosen  as  the  interval  over  which 
the  linearization  and  constancy  of  the  aerodynamic  coefficients  are  assumed 
valid.  The  optimal  control  would  fall  into  the  class  of  dual  control  problems. 
However,  additional  uncertainties  in  the  model  other  than  those  accounted  for, 
unsteady  flow  problems,  as  well  as  the  survivability  dictates  that  the  dual 

I 

control  may  not  be  used.  This  is  because  additional  control  responses  due  to 
the  identification  aspect  of  dual  control  may,  because  of  additional  uncer- 
tainties,  cause  an  extremely  undesirable  response  leading,  perhaps,  to  an 
aggravation  of  the  divergence  problem.  This  may  in  an  extreme  case  lead  to 
the  loss  of  the  aircraft.  Thus,  an  adaptive  open  loop  feedback  controller 
will  be  chosen.  There  are  two  major  techniques  that  may  be  used.  The  first 
in  reference  [7]  yields  the  optimal  open  loop  feedback  controller  for  the 
problem  with  uncertain  parameters.  The  computational  burden  of  the  technique 
is  larger  than  the  second  technique  of  [8]  even  though  it  will  lead  to  better 


perforn'atice  as  it  is  an  optimal  technique.  Since  optimality  certainly  must, 
in  this  application,  include  computational  burden,  the  technique  as  given  in 
1 8]  will  be  extended  to  the  model  following  problem. 

The  partitioning  technique  as  in  [8]  consists  of  using  a control  law 
found  by  solving  the  u-condi tional  control  problem  and  then  weighting  the  p- 
conditional  control  with  the  probability  density  of  u conditioned  on  the 
measurements.  The  technique  includes  the  measurement  conditional  probability 
density  function  for  y in  tne  solution  for  the  control  gains.  This  yields  the 
optimal  open-loop  feedback  control  solution,  but  it  has  the  disadvantage  that 
the  equations  for  the  control  gain  differ  at  each  measurement.  In  this  paper 
a control  law  with  as  little  computational  burden  as  possible  consistent  with 
the  uncertainty  problem  and  good  performance  is  desired.  Consequently,  the 
partitioning  algorithm  will  be  chosen  as  the  adaptive  control  method  in  this 
paper. 


III.  Control  Law  Development 

The  solution  to  the  partitioned  adaptive  control  model  reference  problem 
may  ho  found  by  using  state  augmentation  techniques  The  new  state  c is 
defined  as 


T 

C = 


This  yields  a state  equation  as 

? = A(pj)?  + B(p2)lj 


(6) 


where 


and 


A(ui) 


A(pj ) 
. 0 


B(u2) 


B(p2) 

. 0 - 
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fv 


4 - . 


The  perfonpance  index  n'ay  be  easily  rewritten  as 


= Ejy*  ^|^x^Q(uj)x  + 2u^s(uj,u2)x  + u^R(u2)ujdtj  (7) 


wher  e 


and 


Q(u,)  = [c„a(w,)-a„cJ\[c„a(„p-a^cJ, 


R(u,)  . ♦ Rp. 


(8) 


Thus,  it  may  be  noted  that  the  integral  under  the  performance  index  as 
w(  1 1 as  the  system  dynamics  are  functions  of  and  u^- 

The  partitioned  adaptive  control  law  may  be  found  by  solving  for  the 
ooterininistic  control  gain  conditioned  on  pj  and  ^*1'  '>^2 

^ional  Kdlii’an  filter  estimate,  x(t|  p^  ,P2  » cind  for  the  conditional  density, 

. V^)  where  = ''y(t)  ,t^ £t  j^t) , and  using  as  a control 


u(t) 


-f  f K(tlp^ 

Re  2 


,P2)x(t|  p^  ,P2»'['^)p(Ul 


(9) 


If  p,  and  p-  are  defined  over  discrete  ranges,  then  equation  (9)  may  be  re- 
written as 


^2 


u(t)  = 2 E K(t|pi.,P2.)x(t|pj.,P2.,Y^) 
i=l  j=l  1 J 1 J 


(10) 


where  Pj„((-))  denotes  the  probability 
be  deterr.iin^^d  by  the  solution  of  the 
oroblem,  1 .e. , 


of  the  event  (•)•  The  control  gain  may 
, and  p,  conditional  deterministic 
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K(tlu,  = -R(m,.)'MS{p 


' 1 


i,  ,p„.)  + B(M2.)'p(v^l..M2.,t)]  (n) 


where 


P(w,.,y,.,t)  = P(p,  . ,u„ . ,t)  [A(u  ) - B(u2  .)R(vJ2  .)'■ 

■ S(uj  . ,^2  ■ ) 

- B(u,.)R(u2.)'^S(y.  .,U2.)]'^P(ui.,U2..t) 

- P(pi.»M2.>t)B(u  )R(y  ) ^B(u2  .)^P(vJi  . >^2  . >t) 

H 2j  2j  2j  It 

+ Q(p,.)  - S(mj  . ,y,  .)"^R(y„  .)  ^S(y,  ,^2.) 

v^.  = 1,2 

V.  = 1,2,. ..,12 

J 

with  '"inal  condition 

P(n, . . ,tf)  = 0,  Vi ,j. 

‘1  2j  t 

The  filter  equations  are  the  standard  Kalman  equations 

x(t|)j,,U2.,Y4.)  = A(y..)x(tlvi,.,U2.,Vi) 

+■  B(u  )K(tly,  .,u2.)x(t|y,.,M2.’'^t) 

M 2j  It  2j  c 

+ Kg(tluj . ,P2j)  [y^  ■ Cx(tluj . ,Ji2j  .Y^)] 


wi  th 


where 


x(to|y^.,y2j)  = x(t„) 


Kg(t|yi.)  = V(tly^.)c'^Vp‘l 


(12) 


(13) 


(14) 


w i t-.h 


r 


V(t!u,.)  = A(p, .)V(t|u, .)  + V(t|u  .)A(u, 
M h h M M 


(15) 


- V(tlp^.)C^^"'  CV(tlp^.) 


and 


v(t„!,,.)  = v(t„). 

The  cunditional  probability  density  function  for  y, . and  p,  may  be  computed 
v:a 

Pr(pi.,P2^,Y^)  = 

I *•  t 1 

P(-  ,)P(u  y l|ci<(t|R,..P2..Vt)H^  _,<lt 

Tjl ' "n  ' 

P(  , JP(.  .)exp|y  x(tlPj  .,U2^,Y^)C^Vp‘^y(t^)dty  1 [Cx(t|  Pj . ,p^^  ,Y^)  | | ^ 


, dt ' 


i = 1 , 2 , . . . , 4 1 
j ~ 1»2,...,4.2 

where  P(.-.^)  and  P(j2.)  are  the  a priori  probabilities  for  p^ . and  p^.,  respec- 
tively. The  siirulations  use  the  proper  discrete  form  of  the  equation  as  in 
reference  [9]. 

The  procedure  is>  thus,  to  first  obtain  a set  of  linearized  equations 
where  the  linearization  is  taken  over  the  current  flight  condition.  This 
1 inearization  must  be  updated  frequently  as  the  changing  flight  regime  may 
drastically  cnange  the  dynamic  response.  It  is  assumed  that  over  a given  time 
period  the  aerodvnamic  coefficients  are  constant.  This  assumption  is  valid 


(16) 


over  the  same  region  that  the  linearization  assumption  is  valid.  However,  the 
aerodynamic  coefficients  are  not  exactly  known.  It  is  assumed  that  the  uncer- 
’ain  coefficients  are  defined  over  a discrete  range.  This  procedure  defines 
the  linear  equations  as  in  equation  (1).  The  ideal  model  which  is  dependent 
on  the  angle  of  attack  and  sideslip  angle  is  utilized  along  with  the  measure- 
ment equation  as  in  (2)  along  with  the  determination  of  the  output  equation  (5) 
in  order  to  define  the  remaining  equations  for  dynamic  response.  The  control 
weighting  matrix  R and  model  matching  weighting  matrix  Q must  be  determined. 

r r 

The  final  time,  t^,  must  be  determined.  This  is  the  maximum  time  that  the 
linearization  will  be  assumed  valid. 

Thus,  all  the  equations  necessary  for  control  law  design  are  now  assumed 
to  be  available.  Equations  (13-15)  are  used  to  obtain  the  and  conditional 
state  estimates  for  the  aircraft  state.  This  is  used  in  equation  (16)  to  find 
probability  of  each  Uj  and  ^2-  control  gain  as  in  (11)  is  calculated, 

and  the  control  is  determined  from  equation  (10). 

An  approximate  law  may  be  obtained  by  finding  the  steady  state  gains  K 
and  the  steady  state  Kalman  filter  and  using  these  in  the  control  computations. 


IV.  Ideal  Aircraft  Model 

This  section  uses  the  equations  for  the  A-7  as  given  in  Appendix  A in 
order  to  discuss  the  reasons  for  departure  of  the  A-7  and  in  order  to  yield  an 
ideal  model  with  better  response  in  the  high  angle  of  attack  regime. 

The  equations  of  motion  for  the  actual  A-7  are  given  in  Appedix  A.  They 
are  formulated  using  wind  axes  and  the  aerodynamic  derivatives  are  evaluated 
at  a prestall  flight  condition  of  = 19  deg  and  6^  = 6 deg.  The  incremental 
change  in  rolling  moment  L^.  and  in  yaw  moment  Nl  are  calculated  by 


and 


L'. 

1 


i-ny>xV 


n: 

1 


' - 'bznxb' 


(17) 
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. - A-  • 


r 


1 


vvhere  i deriotcs  trie  particular  state  variable  and  wnere  I and  I are  body 

X y 

axes  nioiixrnli.  of  inertia,  1^^  is  the  respective  cross  product  of  inertia,  and 
L.  and  N-  are  the  aerodynamic  moments  about  the  conventional  aircraft  body 
axes.  As  mentioned  in  reference  HI,  the  major  coupling  which  affects  depar- 
ture for  tne  A-7  is  provided  by  the  kinematic  terms 


and 


^p  = ^0  “o 


^r  = 


Sin  a. 


(18) 


(see  equation  A. ? ) , and  the  aerodynamic  terms  and  The  kinematic  terms 
arise  due  to  the  rotating  coordinate  system.  The  aerodynamic  terms  and 
as  well  as  L'  and  N'  are  given  as  functions  of  a and  B in  Figure  9 of  reference 

b 

[1'.  In  the  reqime  of  a ,,  (the  stall  angle  of  attack)  these  aerodynamic 

S I'd  I t 

ter'ms  change  sign  and,  therefore,  an  aerodynamic  term  stabilizing  at  small  o 
can  contribute  to  the  tendency  for  departure  at  high  a. 

The  A-7  has  a typical  nose  slice  departure.  The  influence  of  the  main 
aerodynamic  derivatives  on  this  departure  is  discussed  at  length  in  reference 
[1|.  Therefore,  in  the  following,  only  points  pertinent  to  finding  the  ideal 
model  are  discussed.  This  discussion  is  based  upon  many  simulation  runs  as 
well  as  physical  insight. 

The  influence  of  the  "effective  derivatives" 


and 


Bo  cos 


Z 

r 


“o 


(19) 


is  not  a major  influence  on  the  ideal  model.  Consequently,  the  ideal  model 
con^ains  these  terms  in  their  original  form.  In  some  of  the  simulations  they 
were  zeroed  and  the  results  were  not  changed  significantly. 

Cor  static  lateral  stability  in  yaw,  N'  should  be  positive.  Figure  9 in 
reference  (11  shows  that  a negative  can  be  expected  for  an  angle  of  attack 
greater  than  17  deg.  'or  a flight  with  high  angle  of  attack  and  small  side- 
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v«n . 


slip,  N'  0 v^in  incr-fcase  the  sideslip  angle  and  will,  therefore,  increase  N' 

a 

w‘’’ch  is,  in  the  high  angle  of  attack,  always  destabilizing.  The  result  is 
tne  high  yaw  rate  of  the  aircraft. 

For  static  lateral  stability  in  roll,  L'  should  be  negative.  This  means 

P 

that  a positive  unwanted  bank  angle  (right  wing  down)  will  induce  a positive 

,i  'esl ip  which  causes  a negative  rolling  moment  with  a decrease  in  bank  angle 

a a result.  A negative  L'  and  a negative  L'  are  the  primary  reasons  for 

p a 

eparture  of  the  A-7  after  a rapid  yaw.  The  is  negative  if  a > 23  deg  and 
ch-  magnitude  increases  with  sideslip. 

■^he  desired  model  for  implicit  model  following  was  obtained  with  several 
[nal . in  iiind.  The  terms  in  the  state  model  were  chosen  with  the  following 
■•esults.  The  yaw  moment  due  to  sideslip  should  be  positive,  > 0.  This 
.1  oecrcase  the  sideslip  and  therefore  the  destabilizing  influence  of  L'  and 
. . .ill?  roll  moment  due  to  sideslip  should  be  negative,  < 0.  An  artificial 
i.t.^tir.  roll  stability  (L^  < 0)  was  introduced  instead  of  the  "natural"  static 
•oil  SMbility  (L'  0)  which  can  contribute  to  departure. 

The  modol  was  chosen  with  the  above  criteria  satisfied.  Simulations  were 


indu;  ■ 

;id  to  help 

( boose 

and  to  verify  the 

model  picked.  These  were 

compa red 

■iCtu;  1 A-7 

in  the 

same  flight  regime 

with  significant 

improvement.  The 

. ' matrix  chosen  for 

the  flight 

condition  about 

“o  = ^9° 

and  = 

6°  is 

-0.063<3  ■ 

-22.68 

0 

-5.766 

0 

0 

3.187 

-32.024 

1 

j 

-0.00G9  ■ 

-0.323 

1 .0 

0 

-0.0995 

-0.0338 

0 

0 

0 

-3.577 

-0.386 

0 

-0.0032 

0.0025 

0 

0 

j 

A = 

in 

0 

0.0122 

0 

-0.1062 

0.3216 

-0.9469 

0.1166 

0.0129 

0 

0 

0 

-1 .0 

-0.849 

0.3323 

-0.5 

0 

0 

0 

0 

1.5 

0.0193 

-0.1276 

0 

0 

0 

0 

0 

0 

1 .0 

0.3397 

0 

-0.0116 

0 

0 

1.0 

n 

0 

0 

0.0104 

0 

ne  available  controls  are  aileron  deflection,  elevator  deflection,  and 
: udOf r deflection.  The  thrust  is  constant  over  the  flight. 

Ihe  next  section  yields  the  results. 
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V . Simulation  Results 

The  cnnfrc  ' lav;  wa'u  used  to  find  tiie  adaptive  implicit  model  following 

Kiv.  ' ■ L'u  A-/  atn,ut  the  « = 19°  and  = 6°  case.  The  model  used  is 

0 0 

d.'i  ned  in  egudt  'un  (L^O).  Measurements  of  all  the  states  corrupted  by  white 
n ise  wern  a.su:iied  available.  The  standard  deviations  for  the  measurement 
■loi  .e  art  as  follows:  velocity  perturbation,  1.25  ft/sec;  angle  of  attack 

perturbation,  0.005  rads;  pitch  rate,  0.01  rad/sec;  sideslip  perturbation, 
0.005  rad;  roll  rate,  0.01  rad/sec;  yaw  rate,  0.01  rad/sec;  roll  angle, 

0.005  rad;  and  pitch  angle  perturbation,  0.005  rad.  In  each  figure  for  this 
case,  the  nomenclature  A-7  corresponds  to  an  open  loop  (6  = *n  = <5,  = 0) 
simulation  of  the  A-7  with  an  initial  yaw  rate  of  -10  deg/sec.  As  is  shown 
in  reference  [Ij  a control  input  typical  of  an  actual  pilot  response  does  not 
control  the  departure.  The  nomenclature  MO  corresponds  to  an  ideal  model 
which  the  control  is  calculated  to  follow,  and  the  nomenclature  MF  corresponds 
the  actual  A-7  response  with  the  closed  loop  control  calculated  in  this 
paper.  Equal  control  weighting  was  used. 

1-  is  assumed  that  the  coefficients  L',  N'  and  L,'  are  uncertain.  These 

u a K 

c./e f f icient'.  nave  a major  impact  on  lateral  directional  stability.  The  true 
value  et  ti»e  parameters  were  3.09,  -1.486,  and  -0.849,  respectively.  It  was 
a .umed  thai  rhf  pos:.ible  parameter  vectors,  {L^^,  N\  Lp]^,  were  1.9,  1.45, 
1.0,  0.55  and  0. I times  the  true  values.  That  is,  the  possible  parameter  set 
was  contained  within  a set  of  five  possible  values.  The  adaptation  took 
place  on  these  three  aerodynamic  coefficients. 

Figures  1-8  show  the  radical  difference  in  response  using  the  control 
law  derived  ir:  tMs  paper.  The  actual  A-7  shows  a buildup  of  roll  rate 
(Figure  5)  followed  by  a rapid  increase  in  bank  angle  (Figure  7).  This  type 
of  behavior  can  indeed  cause  the  loss  of  the  aircraft.  The  responses  using 
the  closed  loop  control  show  that  divergence  is  prevented.  The  response  is 
very  adequate  using  this  control.  Figures  9-11  show  the  control  deflections 
reouired  for  divergence  prevention  for  this  lateral  directional  case. 

Figure  12  shows  the  probabilities  of  each  parameter  being  the  true 
parameter.  It  takes  less  than  1.75  seconds  to  adapt  upon  the  correct  param- 
eter with  probability  0.8. 


Case  two  was  chosen  to  show  the  coupling  between  longitudinal  and  lateral 
dvnan'.ics.  This  case  starts  with  a b deg/sec  pitch  rate  initial  condition 
i^Fiqun.:  13).  Figures  14-16  show  a buildup  in  the  A-7  response  in  the  lateral 
modes  due  to  the  initial  longitudinal  pitch  rate  initial  condition  at  the 
high  angle  of  attack  regime.  Figure  17  vividly  depicts  the  buildup  of  bank 
ingle  as  the  ai>^craft  goes  into  departure  without  the  control  law  used.  These 
figures  also  show  that  with  the  control  law  applied  with  deflections  in 
Figures  18-20  that  the  aircraft  is  prevented  from  departure.  The  control 
laws  in  essence  yield  a soft  decoupling  of  modes  while  controlling  the  air- 
craft. 

Figures  21-28  show  the  standard  deviations  of  the  estimation  error  for 
each  of  the  conditional  Kalman  filters.  The  true  parameter  is  number  3 in 
chese  figures. 

Several  additional  simulations  were  conducted  with  different  flight 
conditions  as  well  as  noise  sequences.  Each  result  is  very  similar  to  these 
typical  results. 


V'  I . Conclusions 

The  control  law  and  philosophy  of  flight  control  developed  within  is 
shown  to  be  an  excellent  method  of  divergence  prevention  in  the  high  angle  of 
attack  regime.  The  control  laws  found  by  finding  steady  state  gains  for  the 
filters  as  well  as  the  control  gains  may  be  readily  implemented  along  with 
the  probability  estimator  for  uncertain  coefficients.  The  control  law  was 
simulated  in  detail  and  shows  excellent  promise  for  control  in  a dangerous 
f I ight  regime. 

The  model  development  philosophy  points  out  many  key  problems  in  the 
high  angle  of  attack  regime.  This  information  of  itself  is  valuable. 
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Case  1 Aircraft  Responses  in  Velocity 
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Fig.  12.  Probabilities  vs  Time 
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Case  2 Aircraft  Resoonses  in  Pitch  Rate 
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Case  2 Aircraft  Responses  in  Roil  Rate 
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APPENDIX 

A 

Equations  of 

Motion 

The  equations  of  motion 

for 

the  aircraft  are 

X = Ax  + 

bu 

(A.l) 

where 

0 

"b 

0 

0 

gpocoso^ 

-ycos  Yq 

'w 

+1 

0 

-Bq  cos 

-BpSinUo 

O 

V.* 

0 

0 

M 

a 

+M-Z 
W a 

Ma 

+M- 

a 

-Po(  ) 

-p“i 

) 

) 

0 

0 

A = 

0 

Y /U- 

a 0 

0 

sin 

-cos 

"‘o 

if-  cos 
0 

if 

0 

0 

L’ 

a 

0 

'p 

L' 

r 

0 

0 

0 

N' 

a 

0 

^B 

% 

N' 

r 

0 

0 
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0 

+1 

tan 

®o 

0 

*^0 

coT^ 

0 

0 

+1 

0 

0 

0 

r 

0 

0 

(A. 2) 

and 
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v/ith  the  state  vector  as  given  in  equation  (1)  and  where  X^- , Y^,  is  the 
k expression  for  )/S^,  is  the  perturbed  total  linear  velocity,  is  the 

free  stream  velocity,  r^  is  nominal  yaw  rate,  p^  is  the  nominal  roll  rate,  g 
. is  gravity,  M.  is  M/n,  is  the  flight  path  angle,  and  6e,  6a,  and  6r  are 

the  elevator,  aile>^on  ana  rudder  deflections,  respectively.  The  aerodynamic 
coefficients  may  he  found  in  references  [11,  12]. 


The  values  for  the  matrices  about  the  a,  6 point  chosen  for  the  simula- 
tions are 


-0. 1025 

0 

0.698 

-0.057 

0 

0 

-2.92 

0 

0 

-0.0037 

0 

0.0255 

-0.292 

0.^31 

1.4 

0.1095 

0.031 

-0.998 

0 

0 

0 

0 

0 

0 

For  the  C matrix  in  the  measurement  equation  (comoare  Eq.  2,  p.  3),  a 
unit  matrix  (8x8)  was  used  in  the  simulation. 


APPENDIX  B 

Model  Following  Program 
-Part  1- 

Control  Calculations  and  Flight  Simulations 

1 . Important  Parameters 

INT ( 1 ) : - not  used 


INT(2): 

(Case) 

1 

open  loop,  ramp  <5^ 

2 

close  loop  0 ->•  6g,  step  0^ 

3 

as  case  2 & u <5j 

4 

open  loop  r(t=0)  = -10  deg/sec 

5 

open  loop  q(t=0)  = 5 deg/sec 

INT(3): 

MODE 

-2 

-1 

run  preparation 

0 

simulation  run 

+ 1 

end  of  simulation 

INT(4): 

NO 

1 

A7,  MO  & MF  simulation 

2 

MO  & MF  simulation 

3 

MF  simulation 

NF 

1 

A7  output  file 

2 

MO  output  file 

3 

MF  output  file 

INT(5): 

0 

A7  derivatives  const 

1 

nonlinear  simulation  (in  older  program 

version  only) 

INT(6): 

KEN 

0 

first  call  of  several  subroutines 

1 

not  first  call  of  several  subroutines 

INT(7): 

I HALF 

max  number  of  bisections  of  the  initial 
ment  in  x - calculation 

time  incre- 

INT(8); 

- 

not  used 

INT(9): 

not  used 
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INT(IO): 

- 

not  used 

R(l) 

XM 

runtime 

R{2) 

DELTIM 

time  interval  for 

new  control 

calculation 

R(3) 

W 

weight  of  control 

R(4) 

integration  time 

for  Riccati 

equation 

R(5) 

time  increment  of 

integration 

(Riccati  equation) 

Ym  = CY  + 
m n 


Subroutine 

CONTROL 


Subroutine 

RKGS 


INT(4)  = INT(4)  + 1 


r 

! Calculation  of  Model 
~ Reference  Control 

L 


^Calculation  of  Y h 
^ with  Subroutine  FCT 
and  Subroutine  OUTP 
• Start  at  DAT(20) 
i Cnd  at  OAT (20)  + DELTIM 

L 


INT(4)  Simulation 

1 A7 

2 


3 


Model 

Model  Fol lowing 


Subroutine 

COTJTjy. 


Subroutine 

SUBCON 


Calculation  of  P with 
I Subroutine  RICCAT 
and  Subroutine  RESULT 


r 

I Calculation  of 
“~l  K 
L 


r 

I Calculation  of 
I V,  X and  KG 
L 
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Subrouti ne 
KALMAN 


Sub>"Outine 

XHAT53 


ENTRY  ^ 


J Calculation  of  V 
’i  and  KG 


I Calculation  of  X 
, with  Subroutine  RKGS, 

I Subi^cutine  XHCAL  (X=^...) 
I^ar.d  Subroutine  XHRES 


53 


r. 

I X 

\ Calculation 

L 


:ine 

B 


r 

.1  I.C.  - Values  for 
I Diagonal  of  V 

•-  i 


I Calculation  of  V \c/ith 

1 Subroutine  VCAL  (V=...) 

I_  and  Subroutine  VRES 


3 . Program  Comments  anu  Source  L_i_s t i n 

The  variable  names  correspona  closely  lo  the  notation  in  the  text.  The 
probability  density  function  (compare  Ecuation  16,  page  8)  was  calculated 
(with  P(p„  ) = 1)  in  the  following  way  (short  notation): 


j=l 


A=/'  xVv,-V„ 


dt  IICXll^  _^dt 

^0  '' 


"T  T -1  -'T  T -1  ' 

X C Vn  Y dt  - / X'C  Vo  CX  dt 
R m / R 


= xVVf^‘^(Y^-  CX)  dt 


To  prevent  overflow  for  longer  simulation  time 


A(i)  = A^(i)  + n(i) 


= const,  A^  < 87.0 


was  introduced.  Therefore, 


^ ” sum(  i ) 


sum(i)  = ^ p|4|  e '' 


and 


A^.  = A^(j)  - A^(i)  + (n(j)  - n(i)) 


The  program  writes  the  results  on  the  disk. 
File  name: 

F0R001.DAT  for  A7 

F0R002.DAT  for  Ideal  Model  (MO) 

F0R003.DAT  for  Model  Following  (MF) 

The  following  pages  show  a typical  dialog. 
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APPENDIX  C 

Model  Following  Program 


-Part  2- 

Graphical  Display  of  Results 

1 . Program  Control 

The  program  was  controlled  by  input  of  a number  (NO): 


NO 

Output 

1 

perturbed  total  velocity 

2 

perturbed  angle  of  attack 

3 

pitch  rate 

4 

sideslip  angle 

5 

rol 1 rate 

6 

yaw  rate 

7 

bank  angle 

8 

pitch  angle 

11 

elevator  deflection 

12 

aileron  deflection 

13 

rudder  deflection 

14 

throttle  (thrust) 

0 

probability  and  standard  deviation 

-1 

as  1 to  8 with  measurements  (noise) 

-2 

end  of  program 
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